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We investigate the mixing of different shapes in the A « 70 region using the adiabatic self- 
consistent collective method in a rotating nuclei. The calculation is done for 68 Se and 72_78 Kr 
which are known to show oblate-prolate shape coexistence at low angular momentum. A pairing- 
plus-quadrupole Hamiltonian, which is known to reproduce reproduce the structure nuclei in this 
mass region, is used. We calculate the collective path between the oblate and prolate minima and 
■ discuss how the collective behaviour of the system changes with increasing angular momentum. 
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I. INTRODUCTION 



The question "what is the shape" of a nucleus is an old one, and probably dates back to the introduction of Bohr's 
liquid drop model Modern nuclear physics has tried to find answers to this question in many areas of the table 
of isotopes. One area of extreme experimental interest is the N « Z region for intermediate mass, where all the 
^\ ' simple models tell us nuclear dynamics rapidly changes with proton and neutron number. Such proton-rich nuclei are 
f**«« . also important in the RP-process of nucleo-synthesis, and are of interest since they are only accessible with modern 
' experimental facilities. 

Mean-field methods are one way to give microscopic meaning to the concept of nuclear shape. They have been 
used quite extensively in nuclear physics @, since they give a good starting point for the description of many but 
the lightest of nuclei. There are a large number of improvements one can make, either at a single minimum in the 
mean- field energy (random phase phase approximation, projection techniques, quasi-particle expansions), or in the 
case that a single minimum is insufficient (time-dependent Hartree-Fock and related methods, Hill- Wheeler method, 
I . etc.). 

In this work we shall apply the theory of large-amplitude collective motion Q, which is a method that tries to 
describe, in a self-consistent fashion, a few-dimensional surface in which most of the low-energy dynamics takes place. 
There are many related methods in use in many areas ofphysics, see e.g. Our method is just one of a large set 
J> | of techniques applied to this problem in nuclear physics [3| . 

Even though it is derived in a very different way, it shows great similarity to the Hill- Wheeler method, which relies 
on the constrained Hartree-Fock Bogoliubov (CHFB) technique. Here a set of Slater determinants determined by a 
small set of judiciously chosen one-body constraints is used to map out what is presumably a dominant part of the 
space of all configurations. The analogy with the technique of generating states for specified angular momentum by 
imposing constraints (cranking) leads to the name "generalised cranking" . The one-body constraints usually consist 
of a few carefully chosen multipole, or particle-hole, operators, as well as a few generalised pairing, or particle-particle, 
ones. For large scale realistic problems such as the description of nuclear fission the number of generalised cranking 
operators needed in order to make a realistic calculation can become very large. There is also no reason to limit the 
constraints to the standard choices; other degrees of freedom, especially those involving spin-orbit interactions might 
also be important. A more satisfactory method should allow the cranking operators to be determined by the nuclear 
collective dynamics itself. 

One such approach, followed in this paper and set out in detail in the review paper Ref. Q as well as in Ref. [5j (a 
similar approach, plus relevant references, can be found in Refs. 

110,1111), leads to a very well-defined approach, 
which can in principle be solved knowing the Hamiltonian and model space. To find the adiabatic collective path 
we use the local harmonic approximation (LHA). It consists of a constrained mean-field problem that needs to be 
solved together with a local random phase approximation (RPA), which determines the constraining operator. This 
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approach lacks practicality, since the size of the RPA problem is, for a system with pairing, proportional to the size of 
the single-particle space squared. Even though enormous matrices can routinely be diagonalised on modern computer 
systems, the algorithm requires repeated diagonalisation of such a matrix, which makes an implementation in realistic 
calculations prohibitively time consuming. A practical solution to this problem was given in Ref. UMM 

It consist 

of solving the RPA in a limited basis of state dependent operators. This was shown to give excellent results for low 
lying RPA solutions in the minimum as well as along the collective path. Our Japanese colleagues have developed an 
interesting alternative based on the full RPA for separable forces, which can be solved by dispersion relations pi El El- 
In this paper we set out to study the the collective motion in rotating nuclei. The method is a gen eralised version 
of the method presented in Ref. where we study collective path at a given angular momentum [lj ■ Experimental 
results show a large variety of shape-coexistence phenomena at low angular momentum in the mass-region around 
72 Kr 0>0>ElEiEl- This is therefore a good testing ground for large amplitude collective motion calculations in 
rotating nuclei. 



II. FORMALISM 

Details of the formalism can be found in our previous work, specifically Refs. [3|, |fj . It is based on time-dependent 
mean field theory, and the fact that this is a form of classical Hamiltonian dynamics. The issue of selecting collective 
coordinates can then be separated in two steps, the determination of the Hamiltonian governing the classical mechanics 
and the approximate decoupling of a few modes in classical mechanics. To facilitate the second task, we usually make 
the additional assumption of slow motion, called the adiabatic approximation. In that limit, it sufficed to expand the 
Hamiltonian to second order in momenta, and we have a problem that can be solved in closed form. The solution, 
as succinctly summarised below, can be stated without any reference to the original nuclear many-body problem and 
the choice of the interaction. The alternative method set out in Ref. @ looks rather different, and relies heavily on 
the use of Slater determinants. As discussed in the appendix, it is very close to our work and only differs in small 
details . 



A. Local harmonic approximation for the collective path in the adiabatic limit 

By parameterising the expectation value of the quantum-Hamiltonian with a set of classical coordinates, which 
corresponds to a classical boson mapping of the density matrix 0], we can parametrise the Hamiltonian locally in 
terms of a set of coordinates and momenta. Clearly, we could also work in the language of differential forms, and 
avoid specifying coordinates, and coordinate frames that are chosen differently at different points 18] . Such elegance 
would hide the technical aspects of obtaining a suitable algorithm. 

We assume a classical Hamiltonian depending on a set of real canonical coordinates, £ Q (a = 1,...,7V), and 
conjugate momenta, np (f3 = 1, . . . , N). The potential V(£) = H(tt = 0) and the mass matrix B a " are given by an 
expansion of H(£,ir) in powers of tt in zeroth and second order, respectively, 

n{^) = v(0 + \B a ^ a -K P + o{^). (1) 

Terms of higher order (such as 7r 4 ) are supposed to be negligible. The central part in our approach to large ampli- 
tude motion is a search for collective (and non-collective) coordinates which are obtained by an invertible point 
transformation of the original coordinates preserving the quadratic truncation of the momentum dependence of 
the Hamiltonian, by 

5" = /"(£), C=9 a (q) (/i,a = l,...,n), (2) 
and the corresponding transformation relations for the momenta p^ and 7r a , 

= 9*^ a , TT a = !%Vv > ( 3 ) 

where we use a standard notation for the derivatives, = -^g a and ff a = -^f 1 . The adiabatic Hamiltonian, 
Eq. QJ, is then transformed into 



H(q,p) = V(q) + -B^p, Pu + G(p 4 ) 



(4) 
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in the new coordinates. The new coordinates g M are to be divided into three categories: the collective coordinate q , 
the zero- mode coordinates q 1 , I = 2, . . . , M + 1, which describe motions that do not change the energy and finally 
the non-collective coordinates q a , a = M + 2, . . . , n. 

The collective coordinate is determined by means of the solution to the local harmonic approach, which consists of 
a set of self-consistent equations. These are: 

1. The force equations 

H, a = Af, a + A I fi a , (5) 

where f 1 are the zero-modes (also called Nambu-Goldstone or spurious modes) and A/ represents a set of 
Lagrange multipliers (which in nuclear physics are sometimes called cranking parameters). A is a Lagrange 
multiplier for the collective mode, stabilising the system away from equilibrium (we shall often call it the 
generalised cranking parameter). 

2. The local RPA equation 

^/^(Ml) 2 /., (6) 

where the covariant derivative V ;a /3 is defined in the usual way (V a /3 = {V a ),i3)- Zero modes correspond to zero 
eigenvalues of the RPA. In principle great care needs to be taken to have zero modes behave correctly away 
from equilibrium. The symplectic RPA |3( is the correct way to do so; unfortunately it is rather cumbersome. 
Below we will discuss an alternative approach to remove the admixture of spurious solution in an approximate 
scheme. 

The collective path is found by solving Eqs. (J5J) and 10 self-consistently, i.e., we look for a path consisting of a 
series of points where the lowest non-spurious eigenvector of the local RPA equations also fulfils the force condition. 
The motion along this line does not necessarily decouple from other degrees of freedom, but we have a "decoupling 
measure" D, that measures how good or bad this coupling is. As discussed in Ref. Q, it measures the change of 
non-collective coordinates for a change in the collective one, which must be zero for exact decoupling. In the minimum 
of the potential the spurious solutions decouple from the other collective and non-collective solutions, and spurious 
admixtures are absent. This is not true for other points along the collective path, especially when approximations are 
made see the further discussion in Sec. Ill CI 



B. Pairing+quadrupole model 

We apply the LHA to the pairing+quadrupole Hamiltonian as described in [2(j • With a constraint on both angular 
momentum and particle numbers the Hamiltonian can be written as 

H' = H - VtN t -luJ, (7) 

r— n,p 

H = E £ 4cfc- E ^{PlPr+PrPl)-^ E QImQ™, (8) 
I r=n,p M=-2 

where ei are spherical single-particle energies, N T are the particle number operators, J is the angular momentum 
operator, Q2M is the dimensionless quadrupole operator 



Q2M = -^Y.( k \ r2Y ^\ l ) 4q. 



(9) 



where 60 = l/V^o is the standard oscillator length and P} is the (dimensionless) pairing operator 

^ = E4d fc • (10) 



fc>0 



In Eqs. CHEJ c\ (c fc ) are the single particle creation (destruction) operators, where fc labels a set of quantum numbers. 
The Hamiltonian J7J is treated in the Hartree-Bogoliubov approximation and it has been shown that at the minimum 
the local RPA for this Hamiltonian is equivalent to the quasi-particle RPA. 
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TABLE I: The interaction strengths used in the calculations. 





k [MeV] 


G n [MeV] 


G p [MeV] 


y8 Se 


0.09894 


0.3437 


0.3441 


72 Kr 


0.091 


0.3469 


0.3514 


74 Kr 


0.08643 


0.3333 


0.3357 


76 Kr 


0.08473 


0.3204 


0.3326 


78 Kr 


0.08452 


0.3077 


0.3287 



The spherical single-particle energies arc calculated from a modified spherical oscillator |2lJ. The interaction 
strengths (see Table []} are chosen to give realistic values of the deformation and pairing gaps in the ground states. 
Our model space consists of two major shells. We follow |20| and multiply each quadrupole matrix element with a 
^-dependent scaling factor. To achieve the same root-mean-square radii for protons and neutrons different harmonic 
oscillator frequencies are adopted for each type of nucleons 20] . The deformation parameters €2 and 7 are calculated 
from the expectation values of the Qq and Q2 operators as 



3 K 
2 Hluq 



£ 2 = ^\/(Qo) 2 + (Q2> 2 , (11) 



C. Removing spurious admixtures in RPA 

Two of the main difficulties of applying the LHA method to realistic nuclear problems are the effort required in 
diagonalising the large-dimensional RPA matrix repeatedly within the double iterative process and the problem of 
spurious admixtures into the collective path. To limit the computational effort we use the method presented in 
Ref. 0, to reduce the size of the RPA matrix. There it was shown that the RPA equation can be solved with 
good accuracy by assuming that the RPA eigenvectors can be described as a linear combination of a small number 
of state-dependent one-body operators. The same set of operators turned our to give good bases for solving the RPA 
problem also at finite angular momentum 12] . To avoid problem with spurious admixture in our collective coordinate 
me remove the spurious admixture from our state dependent basis set. 

We select a small number of state dependent one-body operators F^ k \ k = l,...,n, assuming that the RPA 
eigenvectors can be approximated as linear combinations of the F^ k \ The approximate RPA vector f a is then given 
by 

n 

f, a ^i a = Y, b ^ ) ( 13 ) 
fe=i 

where J 7 ^ is the expectation value of F^ k \ Any spurious admixture to this basis RPA vector is removed by enforcing 
that the vector J-^ is perpendicular to the spurious particle number and angular momentum operators. This new 
basis set is given by 

i 

where are the particle number and angular momentum operators. The coefficients d k can be found by enforcing 
the condition 

£W . oCi) = , (15) 

which leads to a system of linear equations which can be solved for d k . Even though the removal of the spurious 
admixture is of principal importance it turns out to only give a small change to the collective path for the example 
of 72 Kr as studied in Ref. [jj. 



D. Schrodinger equation on the collective path 



After having made a semi-classical approximation, which leads to a classical Hamiltonian, we need to remember 
that we are studying a quantum system. The standard technique to deal with this is to treat the classical Hamiltonian 
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as a quantum one, and to calculate the cigcnfunctions and energies. One must include the zero modes when quantising 
the Hamiltonian, since they describe rotational and other excitations. Quantisation of the Hamiltonian in a metric 
coordinate space turns the kinetic energy into a Laplace-Beltrami operator (see, e.g., Ref. 0) in the relevant space 
and the collective Schrodinger equation can then be written as 

H(X)V(X ) = (X) (16) 

where X represent both the collective coordinate Q and the spurious coordinates. 

Since the potential and the masses are independent of the zero-mode coordinates the wave-function W can be 
separated into various pieces, 

^{Q,<P N ^ P ,n) ^ g-^UiQ)^"^ ^=e^- D{ IK {ny (17) 

y In v Z7r 

where m and k are the quantum numbers for neutron and proton pairing rotation, and /, M, K are the usual rotator 
quantum numbers. We shall be looking at ground states (band-heads) only, and therefore we shall now use I — M — 
K = 0, and since pairing rotational excitation corresponds to a change in particle number, we shall be use m = k = 
as well. Using (|17|l to separate variables, the Schrodinger equation Ijltill can be written as 

-ff- 1/4 ^ (V /2 ^<r 1/4 [/(Q)) + V(Q)U(Q) = EU(Q). (18) 

Since we wish the wave function 'J to be normalisable we require it to be finite, and we must then insist that U(Q) 
goes to zero when g does. In the present work that only occurs when either of the pairing gaps collapses and thus 
B^p N = is zero, and we shall ignore the rotational moments of inertia, which do not change very quickly. Below 
we solve Eq. <|18[) on a grid with the boundary condition that U(Q max ) — U(Q m i n ) — 0. At points where = 
the condition U = holds exactly; for other cases applying this boundary condition will only give an upper limit on 
energy. 



III. RESULTS 



We have performed calculations for a chain of Krypton isotopes in the A = 70 region, to study the effect of changing 
neutron number on the behaviour of these nuclei. To compare with previous work [g, 0J, we have also studied the 
nucleus 68 Se, which will be seen below to be behave in a rather similar way as some of the Krypton nuclei. 



A. 68 Se 



The nucleus Se is of key interest since it has equal numbers of protons and neutrons. It has also been investigated 
using similar methods for J = + by Matsuo et dl. 0,0- Preliminary results of our investigation has been published 
in |23. 

We start by examining the ground state and excited + states in 68 Se where experiment suggests prolate-oblate 
shape coexistence jl-il . We find a collective path going from the oblate minimum over a triaxial energy maximum into 
a prolate secondary minimum, see Fig.^ We see that the main change in the structure along the collective path is a 
almost linear change in the 7-deformation. The energy difference and the potential barrier between the two minima 
is very small, less then 200 keV. We can therefore expect to see a large mixing of different shapes. The collective path 
can be mirrored around prolate and oblate minimum due to the time reversal symmetry. If we were to continue the 
path towards negative values of the collective coordinate Q and towards larger Q the potential would repeat it self 
in a cyclic manner. This collective path is very similar to what was recently found in Ref. 0. The smallness of the 
decoupling measure D gives us reason to believe that we can obtain a very good description of the ground state in 
this way. 

We have solved the collective Schrodinger equation, see Sec. Ill Dl The cyclic nature of the collective path gives 
condition for the radial part of the collective wave function, U, that U(Q) = U(—Q). In Fig. [3] we can see that 
the ground state wave function is spread out almost uniformly along the collective path. The ground state is thus 
neither prolate nor oblate, which should not surprise us since this nucleus is very 7-soft. This results in two (almost) 
degenerate excited states at 2 MeV excitation. This result is similar to those obtained by solving the Schrodinger 
equation on a circle with constant potential. There the ground state is a constant wave function, and the first two 
excited wave functions, cos(f> and sin0, are degenerate, and the energy spacing only depends on the circumference of 
the circle. 
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FIG. 1: Large amplitude collective motion in 68 Se at I = following the lowest RPA solution at equilibrium, a) Energy along 
the collective path, b) The square of the lowest RPA frequencies, c) Deformation £2 and 7/ 180°. d) The dimensionless pairing 
operators (P T )- e) The cranking parameter A. f) The decoupling measure, D. 




FIG. 2: The deformation parameters in the £2—7 plane for Se at I = 0. Labels denote the value of the collective coordinate. 




FIG. 3: The wave function for the large amplitude collective motion in Se at / = 0, following the lowest RPA solution. The 
first three + eigen functions are shown (as red solid, green dashed and blue dashed-dotted line, respectively). These wave 
functions are periodic around Q = and ±1.57. The grey curve (with the scale on the right-hand side) shows the collective 
potential energy for reference. 



7 




Q Q 

FIG. 4: Large amplitude collective motion in 68 Se at / = 2 following the lowest RPA solution at equilibrium, a) Energy along 
the collective path, b) The square of the lowest RPA frequencies, c) Deformation and 7/ 180°. d) The dimensionless pairing 
operators (P T )- e) The cranking parameter A and the rotational frequency u. f) The decoupling measure, D. 



TABLE II: The deformation for the collective states as well as for the mean field at / = and 2 for Se. 



In 


0i 


2 


3 


2i 


2 2 


2 3 


Omfg 


2mf 


£2 

7 


0.211 
25.6 


0.216 
27.5 


0.224 
26.1 


0.250 
48.1 


0.229 
46.5 


0.234 
54.5 


0.254 
60.0 


0.269 
56.4 



At finite angular momentum the behaviour gets somewhat more complicated, as can be seen in Fig. 01 for 1 = 2. 
The potential is no longer symmetric around the minimum and for positive values of the collective coordinate the 
path goes through the triaxial plane into a very shallow prolate minimum. The collective path then continues towards 
negative 7-deformation. The path ends due to an avoided crossing with a higher lying RPA mode, at which point we 
are no longer able to find a stable solution. This avoided crossing indicates a need to include more than one collective 
coordinate. 

This is a generic feature of the method; whenever two modes come closely together, and thus mix strongly, we 
need collective coordinates for both modes to describe collective motion. Fortunately, the decoupling measure gives 
us a clear indication at what point of the collective surface this is the case. States with wave functions that are not 
sensitive to the are of large coupling should be unaffected. 

For negative Q we have a linear increase in 7-deformation until 7 reaches 90° and Q sa —1-7, where there is a 
sudden change in the nature of the collective path as shown by drastic reduction in the neutron pair-field. This is 
accompanied with a slow increase of the 62-deformation with 7 fixed at 90°. Once again the size of the decoupling 
measure tells us that here our single-coordinate approximation breaks down. We can thus not reliable calculate any 
state that has a sizable weight near Q = — 2. 

After solving the collective Schrodinger equation, we find in Fig. that the wave function of the lowest 2 + state 
is concentrated around the oblate minimum with the tail reaching into triaxial and prolate shapes, but with little or 
no weight where D is large. The first excited state has two peaks in the collective wave function corresponding to a 
superposition of prolate and triaxial deformations, and is less reliable due to some small weight where D is large. It 
is only the third excited 2 + state that appears to suffer strongly from weighting where D is appreciable. 

The average deformation for the ground state and low- lying excited states are shown in Table ITTI together with the 
mean- field deformation at / = and 2. We conclude that mean-field and average values can be substantially different 
in these nuclei; also the difference between + and 2 + states is interesting, with the 2 + states much closer to purely 
oblate than the + states. 
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FIG. 5: The deformation parameters in the £2—7 plane for Se at / = 2. Labels denote the value of the collective coordinate. 




FIG. 6: The wave function for the large amplitude collective motion in Se at / = 2, following the lowest RPA solution. The 
first three 2 + eigen functions are shown (as red solid, green dashed and blue dashed-dotted line, respectively). The grey curve 
(with the scale on the right-hand side) shows the collective potential energy for reference. 



If we summarise these results, and compare to the experimental spectra, Fig. [7] we get some interesting results. The 
decoupling measure is small, and up to the lack of fine-tuning of our parameters to individual nuclei, we do expect 
very good results. Thus we predict a pair of almost degenerate + states at about 2 MeV. The results for the 2+ 
states are probably better than we could have expected: the first of these states has little sensitivity to the large 
decoupling measure at Q = —2, but the second has some. The third state must be thought of as totally unreliable. 



B. 72 Kr 

This nucleus has been studied in detail in Ref. 01 • We ag ain start by examining the + states where a prolate-oblate 
shape coexistence has been established experimentally [l4| . We find a collective path going from the oblate minimum 
through a spherical energy saddle point into a prolate secondary minimum, continuing towards larger deformation, 
see Fig. [S] At Q ~ 7.5 there is an avoided crossing between the lowest RPA mode, which we are following, and two 
higher lying modes that are of pairing- vibrational character. This can also bee seen in Fig. [81 as a rapid increase in 
the decoupling measure. For negative Q, at large oblate deformation we see a collapse in neutron pairing, which gives 
a natural end to the path. Nonetheless, even at this end of the path there is an obvious avoided crossing, which 
suggests that the one dimensional approximation is only of limited value for this nucleus. Fortunately, the maximal 
value for D of about 0.6 suggests, based on past experience, that the influence of any second coordinate is small. 
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FIG. 7: Calculated and experimental excitation spectrum for 68 Se as a function of angular momentum. 
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FIG. 8: Large amplitude collective motion in 72 Kr with axial symmetry (following the f3- vibration, the lowest RPA solution at 
equilibrium), a) Energy along the collective path, b) The square of the lowest RPA frequencies, c) Deformation €2 and 7/I8O . 
d) The dimensionless pairing operators (P T ). e) The cranking parameter A. f) The decoupling measure, D. 



In Fig. 13 we plot the collective wave functions and can see that the ground state wave function is concentrated in 
the oblate minimum. The collective wave function has a substantial spread along the collective path and is skewed 
towards the spherical state due to the collapse in the neutron pair-field. Unfortunately, this has substantial weight 
in the area of the avoided crossing. The first excited state has its major component in the prolate minimum with a 
small component in the oblate minimum. The prolate peak in the wave function for the first excited state is broader 
and more symmetric then for the ground state. The third 1 = state is approximate spherical but has substantial 
oblate as well as prolate components. 

At finite rotational frequency the collective path will no longer go through the spherical state, due to the effect 
of the rotational term (it is energetically unfavourable to generate angular momentum in a spherical state). Instead 
the oblate and prolate minim are connected through a path that goes through the triaxial plane. Due to problems 
with level crossings we have to start in both the prolate and oblate minimum, and connect these together to find the 
complete collective path. In Fig. [23 we see that if we start in the almost oblate minimum and follow the collective 
path toward large oblate deformation the situation is very similar as in the non-rotating case above. We first see a 
decrease in the quadrupole moment but after an avoided crossing with a pairing vibration the collective path changes 
its nature and the collective path ends in a neutron pair-field collapse. When we follow the collective path towards 
smaller oblate deformation the system goes through an avoided crossing with a 7-vibrational mode. The collective path 
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FIG. 9: The wave function for the large amplitude collective motion in Kr, following the lowest RPA solution. The first three 
eigen functions are shown (as red solid, green dashed and blue dashed-dotted line, respectively). These wave functions are 
periodic around Q = and ±1.57. The grey curve (with the scale on the right-hand side) shows the collective potential energy 
for reference. 




Q Q 

FIG. 10: Large amplitude collective motion in 72 Kr at 1=2 following the lowest RPA solution at equilibrium, a) Energy along 
the collective path, b) The square of the lowest RPA frequencies, c) Deformation and 7/I8O . d) The dimensionless pairing 
operators (P T )- e) The cranking parameter A and the rotational frequency uo. f) The decoupling measure, D. 




0.1 0.2 0.3 0.4 0.5 4.5 5.5 



FIG. 11: The deformation parameters in the 62—7 plane for Kr at 1=2. Labels denote the value of the collective coordinate. 
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FIG. 12: The wave function for the large amplitude collective motion in 72 Kr at 1=2, following the lowest RPA solution. The 
first three 2 + eigen functions are shown (as red solid, green dashed and blue dashed-dotted line, respectively). These wave 
functions are periodic around Q = and ±1.57. The grey curve (with the scale on the right-hand side) shows the collective 
potential energy for reference. 



TABLE III: The deformation for the collective states as well as for the mean field at 7 = and 2 for 72 Kr. 



in 


0i 


2 


3 2i 


2 2 


2 3 


Omf 


2mf 




0.282 


0.368 


0.150 0.302 


0.403 


0.355 


0.321 


0.334 


7 


60.0 


0.0 


0.0 42.4 


9.8 


19.9 


60.0 


59.4 



will the go into the triaxial plane and follow a path with approximately constant £2 deformation but with decreasing 
7 from just below 60 to just above 0. At this point the collective path goes through another avoided crossing through 
which the computer code cannot find a stable solution. By starting the calculation in the prolate minimum we can 
follow the collective path to the same avoided crossing as as we run into when starting in the oblate minimum. We can 
the add these two path together to get the collective path from the oblate to the prolate minimum and beyond. When 
we solve the collective Schrodinger equation we get a ground state in the oblate minimum which has a substantial tail 
into the prolate minimum as can be seen in Fig. 1121 The first excited state has its main peak in the prolate minimum 
with a second peak in the oblate minimum. 

Unfortunately, for the spectrum of this nucleus, Fig.^l our results suggest that none of the states can be calculated 
reliably, due to the relatively large values of D close to the minimum. Indeed comparison with the experimental data 
is extremely poor; a single collective coordinate is simply insufficient in this case. 

C. 74 Kr 

In 74 Kr the situation is similar to that of 72 Kr. For zero angular momentum the minimum has an oblate shape, as 
can be seen in Fig.^J At large oblate deformations there is a collapse in the proton pair-field proceeded by an avoided 
crossing between a (3- vibration and a pairing vibration that changes the direction of the collective path. Towards 
smaller oblate deformation the collective path goes over a spherical maximum and through two shallow prolate minima 
before it ends in a state with collapsed neutron pairing. The collective wave function, U, of the ground state is shown 




FIG. 13: Calculated and experimental [l4^ excitation spectrum for 72 Kr as a function of angular momentum. 
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FIG. 14: Large amplitude collective motion in 74 Kr at / = following the lowest RPA solution at equilibrium, a) Energy along 
the collective path, b) The square of the lowest RPA frequencies, c) Deformation €2 and 7/ 180°. d) The dimensionless pairing 
operators (P T )- e) The cranking parameter A. f) The decoupling measure, D. 




FIG. 15: The wave function for the large amplitude collective motion in 74 Kr at / = 0, following the lowest RPA solution. 
The first three + eigen functions are shown (as red solid, green dashed and blue dashed-dotted line, respectively). These wave 
functions are periodic around Q — and ±1.57. The grey curve (with the scale on the right-hand side) shows the collective 
potential energy for reference. 



in Fig. E| It shows the wave function is dominated by oblate shapes but contain a substantial tail of spherical and 
prolate shapes. The first excited + state is a mixture of a large prolate component spread out over both the prolate 
minimum and a smaller oblate part. The barrier between the two prolate minima is not enough to build a collective 
state located mainly on each of these. The second excited + state is even more complicated and contain prolate 
spherical and oblate shapes. 

At I = 2 we see in Fig. I16l that the prolate and oblate minimum are again connected via a path that is dominated 
by triaxial shapes. The collective path is ended by collapses in the pair-fields at large oblate and prolate deformation. 
The prolate minimum we see at / = 2 has similar deformation as the second prolate minimum at / = 0. The wave 
function of lowest 2 + state show a strongly collective nature and even though the peak is situated in the oblate 
minimum it has a strong admixture of triaxial and prolate shapes as can be seen in Fig. 1181 The wave function of 
higher-lying 2 + states contain a mixture of prolate, oblate and triaxial shapes. 

In this case, Fig. 1191 we expect very good results for + states since D is small everywhere, and good, but slightly 
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FIG. 16: Large amplitude collective motion in 74 Kr at / = 2 following the lowest RPA solution at equilibrium), a) Energy 
along the collective path, b) The square of the lowest RPA frequencies, c) Deformation £2 and 7/I8O . d) The dimensionless 
pairing operators (P T ). e) The cranking parameter A and the rotational frequency u. f) The decoupling measure, D. 




FIG. 17: The deformation parameters in the £2 — 7 plane for 74 Kr at I = 2. Labels denote the value of the collective coordinate. 



less accurate, results for the 2 + states, since D is slightly larger, 
experimental data. 



Indeed, the results agree remarkably well with 



D. 



Kr 



In 76 Kr the collective path is almost identical to that of 74 Kr. For the unrotated case the lowest minimum has an 
oblate shape as can bee seen in Fig. [3UJ Towards small oblate deformation the collective path goes over a spherical 
maximum and through two shallow prolate minima. At each end of the collective path there is a strong reductions in 
the pair-fields. Due to numerical difficulties of finding a converged solution we are not able continue the calculation 
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FIG. 18: The wave function for the large amplitude collective motion in Kr at / = 2, following the lowest RPA solution. 
The first three 2 + eigen functions are shown (as red solid, green dashed and blue dashed-dotted line, respectively). These wave 
functions are periodic around Q — and ±1.57. The grey curve (with the scale on the right-hand side) shows the collective 
potential energy for reference. 



TABLE IV: The deformation for the collective states as well as for the mean field at I = and 2 for 74 Kr. 



In 


Oi 


2 


3 2i 


2 2 


2 3 


Omf 


2mf 


7 


0.125 
60.0 


0.135 
0.0 


0.098 0.304 
0.0 49.1 


0.320 
43.7 


0.314 
46.3 


0.213 
60.0 


0.259 
57.9 



all the way to the pairing collapsed state. The collective wave function, U, of the ground state is shown in Fig. EH It 
shows the wave function is dominated by oblate shapes but contains a substantial tail of spherical and prolate shapes. 
The first excited + state is a mixture of a large prolate component and a smaller oblate part. The second excited 
+ state is even more complicated and mixes prolate spherical and oblate shapes. 

At / = 2 we see in Fig. 1221 that the prolate and oblate minimum are connected via a path that is dominated by 
triaxial shapes. The lowest minimum is not oblate at all but instead strongly triaxial (7 w 90°). This is due to the 
fact that the collective path goes from the triaxial minimum rapidly through an oblate shaped state to a secondary 
prolate minimum. The collective path ends due to numerical difficulties at relatively small prolate and triaxial 
deformations. The wave function of the lowest 2 + state shows a strongly collective nature and even though the peak 
is situated in the triaxial minimum it has a strong admixture of oblate and prolate shapes as can be seen in Fig. 1181 
The wave functions of higher lying 2 + states contain a mixture of prolate, oblate and triaxial shapes. The energy of 
the states can only be seen as an upper limit since we have not been able to follow the collective path very far in this 
case. 

Thus, see Fig. EH the energies first two + states should be reliable, but the third is somewhat more doubtful, due 
to the peak in D at Q = — 1. The 2 + states are of similar quality. It is thus surprising that theory predicts a low- lying 
second 0+-state that has not been observed in experiment! 




FIG. 19: Calculated and experimental |l7l | excitation spectrum for 74 Kr as a function of angular momentum. 
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FIG. 20: Large amplitude collective motion in 76 Kr at I = following the lowest RPA solution at equilibrium, a) Energy along 
the collective path, b) The square of the lowest RPA frequencies, c) Deformation 62 and 7/I8O . d) The dimensionless pairing 
operators (P T )- e) The cranking parameter A. f) The decoupling measure, D. 




136.00 



FIG. 21: The wave function for the large amplitude collective motion in 76 Kr at / = 0, following the lowest RPA solution. 
The first three + eigen functions are shown (as red solid, green dashed and blue dashed-dotted line, respectively). These wave 
functions are periodic around Q — and ±1.57. The grey curve (with the scale on the right-hand side) shows the collective 
potential energy for reference. 



E. 



For 78 Kr we have behaviour similar to 68 Se in that the lowest RPA mode in the minimum is a 7-vibration. The 
prolate and oblate minimum are therefore connected through triaxial states. The potential barrier between the oblate 
minimum, which is lowest in energy, and the secondary prolate minimum is only 250 keV as can be seen in Fig. 1261 
In Fig. |2H|we see that this leads to a situation very similar to that of 68 Se where the collective ground state wave 



TABLE V: The deformation for the collective states as well as for the mean field at I = and 2 for 76 Kr 



0i 



2 



0, 



2i 



a 



0.048 0.150 0.142 0.253 0.266 0.257 0.211 0.244 
60.0 0.0 0.0 35.3 28.1 34.0 60.0 56.0 
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FIG. 22: Large amplitude collective motion in 76 Kr at / = 2 following the lowest RPA solution at equilibrium, a) Energy along 
the collective path, b) The square of the lowest RPA frequencies, c) Deformation and 7/ 180°. d) The dimensionless pairing 
operators (P T )- e) The cranking parameter A and the rotational frequency uo. f) The decoupling measure, D. 




FIG. 23: The deformation parameters in the £2 — 7 plane for Kr at / = 2. Labels denote the value of the collective coordinate. 



function, U, represent a approximately equal mixture of all possible triaxial shapes with an approximately constant 
£2-deformation and pair-fields. We also see that the excitation spectrum in 78 Kr is very similar to that of 68 Se. We 
have two almost degenerate excited + states that can be interpreted as rotational excitation (in 7). The fact these 
excited states lie at 1 MeV in 78 Kr compared to 2 MeV 68 Se in simply due to the longer 'distance' between the prolate 
and the oblate minimum in 78 Kr. 

At at finite rotational frequency (/ = 2) we have a slightly different picture for 78 Kr. The lowest RPA solution 
in the almost axial prolate minimum is a neutron paring vibration. This leads to a rapid collapse of the neutron 
pair-field at small negative values of the collective coordinate, as can be seen in Fig. [23 At positive values of Q there 
is an increase in the neutron pair-field and a decrease in the 62-deformation. At Q w 1.8 there is an avoided crossing 
that changes the direction of the collective path into the triaxial plane. The path goes towards larger 7 until we reach 
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FIG. 24: The wave function for the large amplitude collective motion in Kr at / = 2, following the lowest RPA solution. 
The first three 2 + eigen functions are shown (as red solid, green dashed and blue dashed-dotted line, respectively). These wave 
functions are periodic around Q — and ±1.57. The grey curve (with the scale on the right-hand side) shows the collective 
potential energy for reference. 




FIG. 25: Calculated and experimental [l^ | excitation spectrum for 76 Kr as a function of angular momentum. 




FIG. 26: Large amplitude collective motion in 78 Kr at I = following the lowest RPA solution at equilibrium, a) Energy along 
the collective path, b) The square of the lowest RPA frequencies, c) Deformation 62 and 7/I8O . d) The dimensionless pairing 
operators (P T )- e) The cranking parameter A. f) The decoupling measure, D. 
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FIG. 27: The deformation parameters in the £2 — 7 plane for Kr at I = 0. Labels denote the value of the collective coordinate. 




FIG. 28: The wave function for the large amplitude collective motion in Kr at / = 0, following the lowest RPA solution. 
The first three + eigen functions are shown (as red solid, green dashed and blue dashed-dotted line, respectively). These wave 
functions are periodic around Q — and ±1.57. The grey curve (with the scale on the right-hand side) shows the collective 
potential energy for reference. 



a almost oblate state with 7 w 60 and Q « 4.1. The collective path passes through this secondary minimum in to 
triaxial plane. At Q w 6.2 we have a second avoided crossing and a transition into a path with decreasing neutron 
pair-field which leads to a neutron paring collapse at Q 7.4. Due to the shallowness of the barrier between the 
prolate and oblate sates the ground state wave function has a substantial tail stretching out of the prolate minimum 
through the triaxial plane into the oblate secondary minimum, see Fig.|M] There are another two low lying 2 + states. 
The first of these is dominated by a oblate component but has also a sizable prolate part. 

In this case, Fig. |23 the spectrum is somewhat similar to 68 Se, due to the extreme 7-soft character of the + 
states. In this case one of the two excited + states has actually been detected. The data obtained in these difficult 
experiments are not of sufficient quality to rule out that this is actually a narrow doublet. The data on 2 + states is 
also compatible with an additional low- lying 2 + that has not been seen in experiment. 



TABLE VI: The deformation for the collective states as well as for the mean field at I = and 2 for 78 Kr. 
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2i 2 2 
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0.397 0.286 
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0.224 
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60.0 


60.0 
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FIG. 29: Large amplitude collective motion in 78 Kr at / = 2 following the lowest RPA solution at equilibrium), a) Energy 
along the collective path, b) The square of the lowest RPA frequencies, c) Deformation £2 and 7/I8O . d) The dimensionless 
pairing operators (P T ). e) The cranking parameter A and the rotational frequency u. f) The decoupling measure, D. 




FIG. 30: The deformation parameters in the £2 — 7 plane for Kr at I = 2. Labels denote the value of the collective coordinate. 



IV. CONCLUSIONS 

We have extended the method of calculating the self-consistent collective path presented in 0, [j| to include con- 
straints on angular momentum and particle number and implemented it for the quadrupole+pairing Hamiltonian 20] . 
The method consists of finding a series of points fulfilling the force equation, where the local direction of the collective 
path is determined in each point by the local normal modes. The local RPA equations and the force equation are 
solved in a double iterative process with constraints on angular momentum and particle numbers along the collective 
path. The method allows us to determine the collective coordinate from the Hamiltonian without having to assume 
a priori which are the relevant degrees of freedom. 
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FIG. 31: The wave function for the large amplitude collective motion in 78 Kr at / = 2, following the lowest RPA solution. 
The first three 2 + eigen functions are shown (as red solid, green dashed and blue dashed-dotted line, respectively). These wave 
functions are periodic around Q — and ±1.57. The grey curve (with the scale on the right-hand side) shows the collective 
potential energy for reference. 




FIG. 32: Calculated and experimental [fH | excitation spectrum for 78 Kr as a function of angular momentum. 



The results confirm the importance of pairing collapsed states for the collective path as suggested in We have 
also seen a very different behaviour of the collective path between the prolate and oblate minimum depending on N 
and Z as well as weather we are looking at the ground state or an rotating state. Without rotation the path either 
goes through a spherical maximum or through the triaxial plane. In the rotating case the two minima are always 
connected via the triaxial plane. This difference comes because at finite angular momentum the system tries to avoid 
the spherical state because the cost of generating angular momentum gets larger as the deformation gets smaller. 
This effect would not have been seen so clearly if the calculation was done at constant rotational frequency instead 
of at constant angular momentum. The difference would be less pronounced if we would have allowed for more the 
one collective coordinate. 

The results shed new light on the behaviour of these nuclei, and make predictions which can be of importance 
experimentally, especially the extreme 7-soft nature of 68 Se and 72 Kr at low angular momentum. Discovery of the 
predicted doublets of excited + states would shed light on the nature of shape mixing in these nuclei. It is extremely 
encouraging that where we expect our method to be reliable (i.e., where D is small) we get a good correspondence 
with experiment, and where we have a poor correspondence, decoupling is poor. We have also seen that the nature 
of shape-mixing and shape-coexistence is indeed quite complicated in the group of nuclei discussed here, and changes 
rapidly along the isotope chain. The picture seems to be less systematic than discussed in the literature. Mixing is 
large in all cases. 

In this paper we have implemented a method to find the adiabatic self-consistent collective path for a nuclei. A 
technique to truncate the basis in which the RPA equations are solved has been improved and a good agreement 
between the full and truncated RPA is found. To solve the RPA equations in a limited basis has proven to be a 
useful and practical way of calculating the collective path within the local harmonic approximation. It remains to be 
investigated how we can include the covariant terms in the RPA in a suitable approximation. 
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APPENDIX A: COMPARISON WITH MATSUO ET AL. 



In the work of Matsuo at ah, the adiabatic large amplitude collective motion (ALACM) is described in what looks 
superficially entirely different than our method. Nonetheless, the final results are so similar that there must be 
similarities. It is therefore useful to analyse the problem in various ways and compare the two methods in a unified 
notation. The key difference is the direct formulation in terms of mean field equations by Matsuo et at, and our 
insistance on classical mechanics. 

Let us initially look at one way to derive the localised harmonic approach in our notation, which we shall than tie 
closely to the work in Ref. |(J. We expand the collective energy (actually the generalised density matrix) to second 
order in small fluctuations, thus obtaining a Hamiltonian. The method is based on "removing the linear terms" , 
which means that at any point P, which is really a time-even HFB state, along the collective path, we approximate 
the Hamiltonian by the quadratic form 

H = \iT a B^{P)^ l3 + V(P) + V a (P)C + \VMP)^ • (Al) 

In order to look at the local fluctuations, we need to remove the term linear in £. We first modify Ti. by constraints 
(such as those for particle number and angular momentum; we shall include particle number only here). With the 
Lagrange multiplier term n(N — Nq) we really have 

H-XN=^it a (B a P - nd na d wp N) (q )^ + V(q ) + (V a - fiN lOt )(q )^ a + X - {V afi (q ) - fiN af3 - M)0 JV i/9 - H,pN a )^ fi . 

In order to remove the linear term, we require that it can be removed by purely collective term, i.e., a term of the 
form XQ, with A defined at the collective point Q = qo- We thus require that 

V a -(MN a = \f ta , (A3) 

with / )Q = Q. a , and we get the quadratic Hamiltonian (where we have absorbed the terms arising from expanding 
(H — fiN) for fixed \x in B and V, and denoted the compound object by a tilde). The Hamiltonian with the constraint 
XQ becomes 

H^^ai^irp + V+^ap-^aN^-^pN^-Xf^fi)^ . (A4) 
It has been argued coherently in Ref. |3| that the final term is equivalent to the covariant term in the derivative 

and we are now left with a simple normal mode problem, which is of course nothing more than the usual RPA. The 
self-consistency is now obtained by requiring f to be one of the normal mode directions. We thus can write the RPA 
as two coupled equation (or as a single one) |23| 

B^f, =Bg a , (A5) 
(V a0 - n ta N, p - Li,pN, a ~ Xf^g 13 = V QQ f a . (A6) 

Let us now compare this approach with that of 0,0, which is cast in the form of a set of variational equations. If 
we take into account that to first order the time-even part of fluctuations is £ Q /\/2, with -K a /\/2 the time-odd part, 
it is relatively straight-forward to show that the three equations (2.34-2.36) in ref. Q would be written as 
(note: 

• we use n where Ref. Q uses A; 

• in our notation A = V(Q),q, where we use the bar to denote a collective quantity, as distinguished from the 
quantity defined in the full space Q . 

• f, a = V25(Q). 
. g> a = V26(P).) 
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V a - vN a - Xf. a = , (A7) 
B^f,p-Bg a = , (A8) 

f, a B af3 l } - Et\ [ f, 7/3 / + f /„] = . (A9) 



The first of these relations is the one that removes the linear term, as above. The second defines the collective 
momentum, and finally the last equation defines the affine connection f, 

Bt= l -Bff. a f^ . (A10) 

Exactly the same relation can be found in Ref. Q . Whereas in the local fluctuation approach we construct fluctuations 
in all directions, the work in Ref. 6] differentiates only with respect to the collective coordinates. This implies an 
assumption about block-diagonality of T that is justified for exact decoupling. We get for Eq. (3.1) in Ref. 0, 

(V a0 - nN^g 13 - V, QQ f. a - V Q [/, q/3 / + f /,«,] - n, Q N a = (All) 

This can also be rewritten as (using Vq = X and /i.q = [i >a g a ) 

(V, a0 - fiN a0 - \f >a p - njN,*)^ = X, Q f, a , (A12) 

which is almost what we would argue constitutes the covariant RPA equation, cf. Eqs. (|A6|I . There is a very subtle 
difference, and that is due to the fact that the matrix on the left-hand-side in lA12l is not symmetric under interchange 
of the labels a and (3. The term that is missing as compared to the local fluctuation approach discussed above 
is —fJL a N t pg^ = —fj, iCe N t Q. The absence comes from treating /i = /j,(Q) in the approach of Ref. [6j as an external 
parameter, whereas, in the language of that paper, we take fi({Q)), i.e., as a function of the mean-field expectation 
value of Q, which clearly changes with the mean field, and thus leads to the additional term. 

We believe that symmetry of all the RPA equations is crucial, and thus are worried about such missing terms. Of 
course, in the current paper we have grossly neglected all covariant terms, which clearly is also objectionable. We 
do wonder, however, if this lack of symmetry is what leads to the if -symmetry non-preserving path seen in Ref. [jj. 
Without angular momentum constraints, the local RPA for axial states will always have eigenvalues with definite 
AK , i.e., to change from following a 3 to a 7 vibration, we must have a real rather than avoided crossing in the RPA. 

The further developments in Ref. |g] are also interesting; they correspond to trying to evaluate f, a p, the variation 
of the collective coordinate operator. Even if it were a one body operator, as we tacitly assume, and is necessary for 
exact decoupling, the equations only fix the 2 quasi-particle part of this operator, i.e., terms going like or cc. 
The approximation made in Ref. [|| is to assume that the c^c matrix elements are identically zero. This is the only 
closed form approximation that can be made, but must be dealt with carefully: it would fail for some problems with 
exact decoupling. In an approach that uses a basis of operators, one could evaluate such terms by the relevant matrix 
elements of the operators, multiplied by their coefficients 

We thus conclude that the underlying expressions in our approach and that set out in Ref. [|| are almost, and 
should be completely, identical. We are concerned about the missing term in the potential expansion. 
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